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ABSTRACT 


Parts I and II of this paper present a comprehensive picture of 
longitudinal wave propagation in a wa r m homogeneous magnetoplasma. 

Part I discusses computed dispersion characteristics for propagation 
perpendicular to the static magnetic field. For a ring electron 
velocity distribution it is found that mode coupling and absolute insta- 
bility can occur. Similar effects are predicted for a spherical shell 
distribution. The Maxwellian distribution gives rise to stable prop- 
agation of undamped waves, and attenuating standing waves. A mixture 
of ring and Maxwellian distributions can give absolute instability 
with stronger growth and lower instability thresholds than for the 
ring distribution alone. Propagation oblique to the static magnetic 
field will be dealt with in Part II [SU-IPR Report No. 326], 
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1 . INTRODUCTION 


In 1958, Bernstein published a comprehensive paper summarizing the 
procedure for obtaining the dispersion relation for longitudinal waves in 
a warm homogeneous magnetoplasma, and describing some of their properties 
Over the last ten years, well over a hundred papers have appeared, treat- 
ing theoretical and experimental aspects of these waves. It is not our 
intention to review this material. That has been done elsewhere (Craw- 
ford 1968). We wish to present systematically in a single source the 
sa lent features emerging from study of the dispersion relation for a 
variety of electron velocity distributions. Many of the results are new. 
Others are scattered through the literature, and can benefit from inclu- 
sion in a coherent picture such as this paper aims to present. We shall 
stress numerical solutions of the very complicated expressions involved, 
-ather than the mathematical formalism, and will interpret properties of 
the predicted propagation characteristics in the light of the refined 
stability criteria that have been developed in recent years (Derfler 1961 
1967, 1969; Briggs 1964). 

The paper divides naturally into two parts. The first deals with 
propagation perpendicular to the magnetic field, for which collisionless 
cyclotron and Landau damping effects do not occur. Part II treats obli- 
que propagation, for which they are important. In each case, similar 
procedure is followed. The dispersion relation for a given velocity dis- 
tribution is first solved for complex frequency and real wavenumber, to 
establish whether instability is predicted. If so, stability criteria 
are applied to determine whether the instability is convective or abso- 
lute (Sturrock 1958) . Finally, if absolute instability is absent, the 
dispersion relation is solved for real frequency and complex wavenumber. 

The plan of Part I is as follows. Section 2 describes first the 
procedure for obtaining a dispersion relation, then comments on the con- 
ditions under which instabilities may be encountered, and finally pre- 
sents the stability criteria used to classify them. Sections 3 to 6 ana- 
lyze the propagation for specific electron velocity distributions. Sec- 
tion 7 provides a brief discussion of the results. 



2 . THEORY 


The basic equations from which warm magnetoplasma propagation char- 
acteristics are derived in the quasistatic approximation are Poisson's 
equation and the linearized Vlasov equation (Stix 1962), 

P_( £jt) n e - 

V . E(r t) = — — / dv f (v r t) , (l) 

0 0 J 
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where the plasma is assumed uniform and infinite, with a time invariant 

electron velocity distribution, f.(v ,v ); v and v are the compon- 

’ ° J-’ II 1 || 

ents of the electron velocity, v perpendicular and parallel to the 
external magnetic field, B^; f^ is the lowest order term in the pertur 
bation expansion of the electron velocity distribution; E is the pertur 
bation electric field; is the external charge density representing 

the source of the perturbations; -e is the electron charge; m is the 
electron mass; n^( = n e Q = h^q) * s avera 6 e charged particle den- 
sity, and is the permittivity of free space. Ion motions will be 

neglected . 

To solve (l) and (2), we introduce a Fourier transform in space, and 
a Laplace transform in time, defined with their inverse transformation as 


E(k,o>) = / dt / dr exp i( k • r-to t )e( r , t ) 


(3) 
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where C is the Laplace contour in the lower half complex CD-plane. 
Transforming (l) and (2), and combining, yields 


E(k,CD) 


e o AKlO 


(5) 
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where k = + k and the plasma equivalent permittivity is, 



where k^ and k^ are the components of the wavevector, k, perpendi- 
cular and parallel to the applied magnetic field; <A3^ = n^e^/e^m)^ ' ls 
the electron plasma frequency; £U c (= eB^/m) is the electron cyclotron 
frequency; is the nth order Bessel function of the first kind, and 

03^ is the imaginary part of u>(= + icu^). 

Equation ( 5 ) can be inverted with the aid of (4) to yield, 


E (r, t ) = f§ ex P( i£Ut )£(“ .£) > 


G(o), r ) = / 


dk exp(-ik*r)o (kffl)k 


( 2 ^) : 




( 7 ) 


Dispersion Relation 

The dispersion relation is simply K(u3 > k) = 0, and describes the 
propagation of waves varying as exp i^t-dt^r). If we solve for real), 

exponential growth and decay are implied for 03^ < 0 and U3^ > 0, respec- 
tively, while 03^ = 0 indicates undamped propagation. We shall return to 
this below, but will first specialize (6) to the case of interest in 
Part I, i.e. perpendicular propagation (k^ = 0). We obtain 
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For computational purposes, an integral representation is useful, 
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where Q has been written for cd/cd Note that (8) and (10) are analy- 
tic over the entire complex CD-plane, except possibly at CD = no^. 

Although detailed numerical studies are generally needed to appre- 
ciate the implications of (10), consideration of some simple limiting 

2 2 

forms may be helpful. For example, consider the case where cd^/cd^ « 

and 0 ) « ncD . We may approximate the infinite series in (8) by the nth 
c 

term, so that 


or 


ncr 


K(ci) k ) «* 1 - a (k ) — = 0 . 

* x ...2 n' i (D — ncD 


P n' x ' cd — ncD 

CD C 
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Solving for the frequency, we find, 


cd( k ) ^ ncD 
v x c 


CD 

1 + “f 'n<\> 

0 ) 
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( 12 ) 


( 13 ) 


establishing that there is a mode near each harmonic of the electron 
cyclotron frequency. It may be shown more generally that there are always 
cutoffs (jt = 0) at ncD c (|n| > l), and the upper hybrid frequency 
(=(cd^ + a^)£), and resonances (jc = °°) at na) c (|n| > 0). For this 
reason, the waves are often referred to as 'cyclotron harmonic waves' 

(CHW) . We shall adopt this title in what follows. 

Stability Considerations 

Two topological forms of a n ( k i ) will be examined in Sections 3 - 6 : 

first, such that a (k ) > 0 for all rea] k , and second, such that it 
* n x 1 


h 



undulates about the line a n ( k i ) = 0- Baldwin and Rowlands ( 1966 ) have 
shown that the first case is a sufficient condition for stable CHW propa- 
gation. It is satisfied by distributions such as the Maxwellian with 
df^/dv^ <0 (v > 0). It follows that a necessary condition for insta- 
bility is that Sf^/dv^ > C for some range of v^. 

It will be noted from(ljj) that if a (k ) can undulate, the ampli- 

2° 2 A 

tude of the undulation increases with 03^/03^. Consequently, the loops 
above a given harmonic may ultimately intersect those below the next 
harmonic. This mode coupling will be shown in Section 3 to produce 
instability, but it may be appreciated by making a two-term approximation 
to (8), 


03 
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This can be rearranged as the quadratic equation, 
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03 - q 03 03 + q 03 = 0 . 

1 c 2 c 1 


03 _ 

- (2n + 1) [„. n + (n 4- 1) %+I j , 


q 0 = n(n + l) 


03 

1 + — £ (a + a ) 

03 2 11 0+1 
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(H) 


(15) 


Complex solutions, and hence instability, can occur if there are values 

p 

of k for which q^ < l4q.,. If this turns out to be the case, either 
within the validity of ;he approximation leading to (l5), ° r after study 
of the full dispersion relation, a more detailed examination is required 
to determine the nature of the instability. 

Instability Classification 

CHW are solutions of our basic linearized differential equations, 
varying as exp i^t-Jk.^r) , with 03 and £ related via the dispersion 
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relation K(co k) = 0. If k is real and to. < 0, the wave amplitude 
will increase exponentially with time for any r, no steady state condi- 
tions appear possible. It was pointed out by Sturrock (l9i;8), however, 
that there are two distinct types of instabilities, one of which does 
lead to a steady state. The difference becomes apparent only after a spec- 
trum of waves is superimposed by carrying out the inversion described by 
(7). As t -» 00 , ( 7 ) may take on one of two forms when instabilities are 

present: the amplitude of the fluctuating electric field may grow in 

time at every point in space, or the growing disturbance may propagate 
away, leaving the plasma in a quiescent condition. The former case is 
referred to as 'absolute' instability, and the latter as 'convective' 
instability. The importance of making this distinction should now be 
clear. If the plasma is convectively unstable, a steady state exists, 
and it should be possible to excite sinusoidal oscillations, characteri- 
zed by 00 real and k complex. 

The nature of the instability is determined by the analytic proper- 
ties of G(oo r) in the lower half complex 00 -plane: analyticity there 
permits convective instability, while a singularity implies absolute 
instability (Derfler 1961 , 1967, 1968; Briggs 1964). The singularities 
occur at frequencies where the contour of integration (or surface of 
integration if the integral is multidimensional) of G(oo,r^) is pinched 
by zeros of K(oo,k). For example, for perpendicular propagation as 
exp i(oot-k^x) we can write, 


f dk exp( “ik x)p (ct),k ) 

£.<“>£> s £ G (“’ x) = £j — e n k ) > 

** 0 1 , 1 


( 16 ) 


where ic is a unit vector along the x-axis. Figure 1 shows the contour 
of integration that uniquely defines G(co,x) over the entire complex 
co-plane. As to moves away from the Laplace contour, it may be necessary 
to deform the contour into the complex plane to avoid an approaching 
singularity. This deformation becomes impossible if zeros of K(cn,k^) 
pinch the contour as illustrated. 

It can be readily demonstrated that at the print of pinching 
dK/dk^ = 0, or equivalently *co/ak^ = 0, implying that G(co, x ) has a 
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branch point at the same frequency as k^(co) has. We may make use of 
this fact to provide an alternative to the pinching-pole criterion just 
described. If a mapping of the real k^ axis into the complex co-plane 
traces a loop located entirely in the lower half of the plane, a branch 
point of k^(oo) must be encircled. This implies that G(co,x) has a 
branch point at the same frequency, and hence that the plasma is abso- 
lutely unstable. We shall use this later to show that all of the insta- 
bilities described in Part I are absolute. 
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RING DISTRIBUTION 


The ring distribution describes monoenergetic electrons, moving 
only in the plane perpendicular to the magnetic field, and uniformly 


distributed in velocity space on a circle with radius 
written in terms of Dirac 5-functions as 


Oi' 


It may be 


f (v v ) = ~ — - — 5(v - v n ) 6(v ) 

O' i» || 2 k v Q ^ l Oj/ || 


With this distribution, integration of ( 9 ) and (ll) yields 



f o ( t) = J o(^x sin i) * 


(17) 


(18) 


where p^ = k^v . Substitution of (l8) in (8) and (ll), and some 

manipulation, gives the following alternative forms for the dispersion 
relation, 


K(cn, ki ) 


= 1 


^ ^ .. dxT^p ) ncD 

1 n Vh i' c 
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...2 L , u 
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c n=-oo ll c 


2 it 


CD r 

, p / , sin Qt ■, /^ T \ 

- 1 + ^ J dT nirFh sin T J o (^x cos 2 ) > 


CD 


1 - i ^ J >x>' 


= 0 


CD 


(19) 


Dispersion Characteristics 

Figure 2 shows numerical solutions to ( 19 ) for several values of 

00 2 2 

CD /cd . We note that for cd /cd < 6.62 the branches undulate about suc- 
p c pc 

cessive cyclotron harmonics. It is easy to see from the final form of 
(19) that if Q = n, then J n (p L ^)a j n (pi^) /Su. ^ = 0, implying that the 
modes cut the cyclotron harmonic lines at the zeros of J n (pc ^ ) and its 
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S3 2 

Fig. 2 


Vo± /w c 


. CONTINUED. 


11 


derivative. 


As CD^/CD^ increases, the loops above a given harmonic approach the 

loops below the harmonic immediately following it. The points at which 

the loops can couple must always lie between ct and Co , \ . where 

nm (n+ljm’ 

Q represents the mth zero of J (u ). The values of od-/ 0)^ at which 
nm n' r i p c 

the first few modes just touch are given in Table 1 . 


Distribution 


Frequency 

Band 



OCCD/CD <1 
c 

1< CD/CD <2 
c 

2<cd/cjd <3 
c 

3 < cd/cd < 4 
c 

Ring 

17.02 

6.81 

6.62 

6 . 9 h 

Spherical Shell 

— 

215.38 

57.05 

U 7.91 


2 . 2 

TABLE 1 . Instability threshold values of CD /cd 
P_c 


After coupling has occurred, there are ranges of p^ in which 

purely real solutions for cd do not exist. The real parts cd^/cd^ 

the complex conjugate roots for frequency are indicated in figure 2 by 

fine lines, and the corresponding imaginary parts CD Vcd^ are shown 

dotted. Note that these growth rates can become very strong indeed. 

For example, when cd /CD ~ 1 growth rates of the order of 50 dB per 
’ ic 7 

cyclotron period (= are implied. 

An interesting feature of figure 2 is the occurrence of an insta- 
bility with CD /CD = 0 , implving that plasma fluctuations will grow in 
r c ’ 

time without propagation. The threshold conditions for this are obtained 
by setting CD = 0 in ( 19 ), and expressing the result in the form 



J oW> 


( 20 ) 


Since p is real, this equation can be satisfied only if the line 
-1/2(cd‘"/cd^) lies above a minimum of J^(p^)/n^. When this occurs, 

there exist purely imaginary roots of the dispersion relation. This can 
occur only if p lies between succeeding pairs of zeros of and . 
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Stability Analysis 


Since complex^ solutions occur with U3 < 0 for k real, we must 

i x ’ 

apply the instability criteria summarized in Section 2. Figure indi- 
cates the use of the pinching pole criterion for the particular case of 

cn"' /o)^ = 20 and 1 < 0) /U3 <2. When contours A, B, C and D ore mapped 

pc r c ’ ’ 

into the complex K^-plane via (19), a saddle point is found at 

= ( k 1 v Q i ^ u c )o = + * 0.18. The corresponding frequency is 

(co/03 c ) o = 1.29 - i 0.47 . Since the saddle point is formed by the merg- 
ing of two roots of the dispersion relation that originate on opposite 
sides of the real axis, and hence the contour of integration, pinching 
occurs, and G(cd^x) is singular in the lower half complex U3-plane at 
U3^. Hence, the plasma is absolutely unstable. 

To illustrate the loop criterion, figure Jb shows the type of locus 
followed in figure Ja by the complex 03 roots of the CHW dispersion rela- 
tion as k increases along the real axis (k . = 0). All the instabili- 
l Xi 

ties covered in Part I are of the mode coupling type, and exhibit similar 
behavior. To separate the two branches in figure Jb, a second locus has 
been drawn along which k = 6, where 8 is small and positive. The 
equation describing this locus is obtainable by expanding the frequency 
03( k^ r + i5), in a Taylor series, 

Su3( k ) 

U3(k 1- i&) %03(k ) + iB , (21) 

v xr xr Ak ’ 

xr 

where 03( k ) and the sign of its derivative are evident from dispersion 
' xr 

diagrams such as figure 2. A loop in the lower half U3-piane is invaria- 
bly indicated. This leads to the conclusion that when CHW instabilities 
occur for perpendicular propagation, they are absolute. 
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(a) Pinching pole criterion 



(b) Frequency loop criterion 


Fig. 3. APPLICATION OF INSTABILITY CRITERIA TO PERPENDICULARLY 
PROPAGATING CYCLOTRON HARMONIC WAVES. 
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SPHERICAL SHELL DISTRIBUTION 


Even if a pure ring distribution could be set up instantaneously, 
it would diffuse in velocity space due to collisions, and nonlinear wave/ 
particle and wave/wave interactions. In this section, we shall consider 
as a relevant standard form the spherical shell distribution, 




) = 


'Jtv, 


“ 5 < v - V 


( 22 ) 


v = (v 2 + v^)^. 
1 II 


where 
have the same speed, 


This describes a plasma in which all electrons 


0' 


and are isotropically distributed in velocity 


space on a sphere of radius Vq. The ring distribution would, in princi- 
ple, relax to this case after the electrons have undergone collisions in 
which only momentum is transferred. If energy is also transferred during 
the collisions, the distribution will evolve toward the Maxwellian, which 
we shall study in Section 3 . 

With (22), integration of ( 9 ) and (ll) yields 


a (k ) = 
n' i' 


5 2 " 


F 0 (t) = 


sin 2£ sin — 

2 

2? ^sin | 


(23) 


where 5 = k v_/<d . Substitution of ( 23 ) in ( 8 ) and (ll), and use of 

= 1 i U c 

some Bessel function identities gives the following alternative forms 
for the dispersion relation. 


K(cu,k x ) = 1 


CD 

- 4 V 




ncc 


CD 


L-, 

n=-oo 


CD - nCD 


= 1 


LU - 

+ i/ 


~ sin Qt sin — sin^2§^ cos — j 

dT = ^ — 


5 sin Qn 


2 r 


CD 

= 1 + \ “2 
e cd 

• 1 C 




= 0 


( 24 ) 


where the final form is expressed in terms of Lommel’s function. 
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n/2 

S -l,2t/ 2 V * ,-0 7lnf to / dT cosalT cos ( 2 ? 1 cobt ) • (25) 

0 


Dispersion Characteristics 

2 ^ 

Figure 4 shows numerical solutions of (24) for values of CL'~/U}" 

pc 

similar to those used for figure 2. The modes cut the cyclotron har- 
monic lines at the zeros of J (2? ). This may be seen most easily from 

2n 

(25), in which the integral reduces to jtJ n (2£^)/2 at = n. The 
undulations are less violent, and the instability thresholds may be com- 
pared in Table 1 with those for the ring distribution. They are much 
higher, and we note that no zero-frequency instability occurs. This 
last point may be proved as follows. 

To determine the threshold for instability, we put to = 0. The first, 
form of (24) then reduces to 


1 + 


0)^’ 

_E 

c 2 


1 - 

f 2 

1 


= 0 


( 26 ) 


Since the left-hand side of this expression is positive definite for real 
and nonzero £ , instability can never set in. 

When there is mode coupling and instability in the higher passbands, 
a similar analysis to that of Section 3 shows the instability to be abso- 
lute. The generally weaker growth rates are attributable to the spread 
in transverse energies associated with the spherical shell distribution. 
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Fig. 4. DISPERSION CHARACTERISTICS OF PERPENDICULARLY PROPAGATING 
CYCLOTRON HARMONIC WAVES: SPHERICAL SHELL DISTRIBUTION. 
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5. MAXWELLIAN DISTRIBUTION 


The final basic distribution to be considered is the Maxwellian, 


f o ('/>[- 


ic:jt v 


V2 / v c + v c 

exp [- — — 11 , , 

2v 


( 27 ) 


wher». v^= ( xT^/m^) 2 \ is the electron thermal velocity. The now famil- 
iar procedure gives 


exp( -X)l (X) 

» (k ) = 

n' l X 


, F q (t ) = exp(-2X sin 2 tt ) , 


( 28 ) 


where X = (k^v^/c^) • The dispersion relation may be expressed in the 
alternative forms. 


K(®,k A ) 


= 1 


^ «— > exp(-X)l (X) no) 
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U J w y 

-i 1 

OJ n — — 0 


n=-oo 


co - nca ’ 
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2 it 


<o - 

, p / , sin Or / 2 t \ 

= 1 + —T I dr — — — sirrr exp -2X cos — ) = 0 . 

^2 J sin flit \ 2 / 


(29) 


Dispersion Characteristics 

Figure 5 shows numerical solutions co( k real) of (29) for several 

2 2 ^ 

values of co^/co” 1 . Since a (k ) > 0 for all k real, no mode coupling 
pc n' i' l ’ 

occurs, and the propagation is stable. Experimental observations of CHW 
propagation show that these theoretical dispersion characteristics are 
approached closely in the laboratory (Crawford et al 1967). 

Since a Maxwellian plasma is free of absolute instabilities, steady 
state solutions of ( 29 ) can exist for CD real, and most generally k^ com- 
plex. Figure 6 illustrates this: In addition to the ci real, k^real 
branches already shown in figure 4, there are the C) real, k^ complex 
branches numbered 1-6. In reality there are an infinite number of them, 
and a detailed study has been made by Buckley ( 1968 ). In general, their 
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k l v t /u c 


Fig. 6. DISPERSION CHARACTERISTICS OF PERPENDICULARLY PROPAGATING 
CYCLOTRON HARMONIC WAVES: MAXWELLIAN DISTRIBUTION (k^ complex). 
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spatial decay rates are so high as to make them uninteresting. 

These solutions of (29) pose a minor paradox, since they suggest 
that damped waves can propagate in a lossless plasma. Similar situations 
occur for other modes of propagation, e.g. surface waves on a bounded 
plasma column (Ristic et al, 1969), so it is worth pointing out the cor- 
rect interpretation of this curious result. 

Consider (l6). If the integration is carried out, defining the con- 
tour as in figure 1, we obtain the following representation of G(cd,x), 

N V ex P(- ik n x)o sK k n 

CU,X ' 2^ ie_k BK(to k )/ak 

0 n ’ n x 

where k (to) is a root of k(oj k ) = 0, and the summation extends over 
n ’ n ’ 

all k that have a negative imaginary part when to is on the Laplace 
n 

contour. When cu is real, the terms in this series with complex k 

’ n 

can be put in pairs, since the symmetry of K(co,k ) implies that if k 

* 

is a root with a negative imaginary part, then so is -k , where 

7 n 7 

denotes complex conjugate. Thus, the series is constructed of paired 
expressions 



exp( -ik x)p k ) 
v n s ’ n 

ie k ?lK(ffl k )/Bk 
0 n ' * n i 


exp(ik*x)p (u>, -k*) 
v n s v ’ n 

ieJk*dK(to -k*)/dk 
0 n ' ’ n x 


exp( -ik x)p (<u k ) f"exp(-ik x)p (a> k ) 
v n p s v ’ n ' n /h s ’ n 

= ie^k ^K(o) k )/3k ie^k SK(o) k )/3k 

0 n ' > n' x 0 n ’ n' x 


= - 2G (co k ) exp( k .x) sin(k x-6) , (jl) 

nO , n ' ni nr ’ 


where G ^(cuk ) and 6 are defined by, 
nO ’ n ’ 


p (^,k n ) 

“no' 10 -^) exp(i9) * e k a K (m k )/ak 

0 n v * n n 


(32) 


and the source function, p (<h,k ), is real and symmetric in k . Since 

S X X 

x > 0 and k . < 0, the implication of ( 31 ) is that the U) real, k com- 
ni 7 1 

plex solutions of ( 29 ) are associated with the excitation of standing 
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waves that decay exponentially away from the source. Damped propagating 
waves will not be excited. 


Collisional Effects 

For isotropic distributions such as the Maxwellian, momentum trans- 
fer collisions may be taken into account approximately by suitably modi- 
fying the dispersion relation. To do so (Tataronis and Crawford, 1965 ), 

p P 

first replace cu by (tn-jv), and then by 1 - jv/to), where v 

is the effective momentum transfer collision frequency and is assumed 
independent of speed. With these substitutions in (29), the four most 
lightly damped curves of figure 6 become as shown in figure 7- Two qua- 
drants have been drawn for k v /ou . to demonstrate clearly how the 

lr t c ’ 

real parts of the curves join. The imaginary parts indicate strongest 
damping where the group velocity (dtn/dk) of the real branches is low 
or zero for v = 0 . It should be noted that the attenuation is increas- 
ingly heavy in the higher passbands, and may easily be of the order of 
tens of dB/wavelength for the value of chosen. 
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6. MIXED DISTRIBUTION 


It is frequently the case experimentally that a group of energetic 
electrons interacts with a stable background plasma. Since the ring and 
Maxwellian distributions represent extremes of instability and quiescence, 
it is natural to investigate a situation in which proportions of these 
two distributions are mixed. We shall do so in this section with the 
aims of determining the parameter ranges over which instability will 
occur, and the strengths and nature of these instabilities. We combine 
( 19 ) and ( 29 ) to obtain 


K(®,k A ) 


= 1 


CD 2 

_ _E 

cd 2 


Q 


^ exp( -X ) I n ( X ) 


ncD 


CD - nCD 


n =- 00 


v 1 

+ (l-Q > J--B -t- 


new 




CD 


nCD 


n =-t» 


= 0 


( 33 ) 


where a defines the proportion of the total electron density that is 
Maxwellian. 


Dispersion Characteristics 


We shall not attempt a detailed review of dispersion characteristics 

2 2 

in the CD - k form, since we now have not only cd /cd but also a and 
x » P c 

v /v as variable parameters. We know the limiting forms for Q = 0 
Ox t 

and 1 (figures 2 and 5), and can deduce that reducing a from unity will 
cause undulations, due to the ring component (l-a), to appear on the basic 
Maxwellian curves. The "wavelength'' of these undulations along the k^ 
axis will depend on v^ /v . We will return to a specific example shortly, 

°x t 2 2 

but will investigate first the threshold values of cd /cd at which mode 

P ^ 

coupling due to the undulations will occur. Computations based on (33) 
yield the data shown in figure 8 for the first four passbands. 

When Q = 0 (ring group alone), we retrieve the instability thres- 
hold values of table 1. As a increases, the threshold conditions 

2/2 

become strongly dependent on v /v . In the first passband, 


2k 
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Fig. 8. PERPENDICULARLY PROPAGATING CYCLOTRON HARMONIC WAVES: 
INSTABILITY THRESHOLD FOR A MIXTURE OF (a) MAXWELL IAN AND 
(1-a) RING ELECTRON VELOCITY DISTRIBUTIONS. 
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Fig. 9. PERPENDICULARLY PROPAGATING CYCLOTRON HARMONIC WAVES: 
MAXIMUM INSTABILITY GROWTH RATES FOR A MIXTURE OF (a) MAXWELLIAN 
AND (1-a) RING ELECTRON VELOCITY DISTRIBUTIONS. 
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increases monotonically with v^ /v t> The remaining three bands have a 
more complicated structure, however. For a given a, note that the 

instability threshold may be lowered if v. /v is sufficiently large, 

Ul 1 2 2 

i.e. instabilities are predicted for values of cu'/U) that correspond 

pc 

to stability when the electrons are exclusively in the ring group. A 
concomitant of this is that the maximum growth rate associated with the 
ring distribution alone is increased by addition of some Maxwellian 
plasma. This is illustrated for the passband 1 ^ s 2 in figure 9. 

Study of figure 8 indicates that for v^ /v^ “* 00 there is instabil- 
ity, even for a vanishingly small population in the ring distribution 
(a l). The threshold values corresponding to this limit are 

= 5 (l <: oj/to S2), 8 (2 s; <; 3), 15 (3 <; cd/uj^ <; 4), which 

we recognize as the values for which the upper hybrid frequency equals 

. Near these frequencies propagating branches of the Maxwellian 
c 

plasma can easily be coupled together by the undulations caused by the 

presence of some electrons in a ring distribution. This point is made 

2.2 

explicit by figure 10 for the case 0 )^/ 0 )^ = 3 , v Qj/ v t = ^ ’ which is 
relevant to the diagram for 1 £ ^/ £J0 C ^2 in figure 8. As Q decreases 
from O.96 to 0 . 90 , we spp from figure 10 that mode coupling takes place, 
and there is absolute instability. The growth rate will maximize with 
decreasing a, then decrease again. The final sequence shows how the 
modes finally become uncoupled as a decreases from Q = 0.40 to O.3O. 
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7. DISCUSSION 


In Part I of the paper, we have explored the dispersion character- 
istics of perpendicularly propagating; CHW for u set of basic electron 
velocity distributions, and it is convenient to summarize here the 
principal results. The salient features of the ring distribution are 
its undulatory propagation branches, and the mode coupling leading to 
very strong absolute instability. The spherical shell distribution 
shows weaker versions of the same phenomena, but is not subject to a 
zero- frequency instability. The Maxwellian distribution is associated 
with stable propagation, and leads to either purely real propagation, 
or attenuating standing waves. There is no collisionless damping. 
Finally, our study of a mixture of ring and Maxwellian distributions 
indicates that the presence of a Maxwellian background may cause abso- 
lute instability with higher growth rates, and at lower thresholds 
2 2 

(oj fjj ) than for the ring distribution alone, 
p c 

We may speculate on the results to be obtained for other basic 
velocity distributions and mixtures. If the shape is given in analytic 
form, it may be possible to determine immediately from the form of 
a^(k ) whether there should be instability or not. If so, its strength 
might be inferred crudely by considering the transverse velocity spread 
relative to those of the distributions studied here. It would be 
academic to push this too far, however, as we shall show in Part II 
that the instability picture becomes much more complicated when oblique 
propagation is considered. In particular, growth rates may be higher 
and instability thresholds lower. 
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